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The Newtonian limit at intermediate energies 
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We study the metric solutions for the gravitational equations in Modified Gravity Models (MGMs) . 
In models with negative powers of the scalar curvature, we show that the Newtonian Limit (NL) is 
well defined as a limit at intermediate energies, in contrast with the usual low energy interpretation. 
Indeed, we show that the gravitational interaction is modified at low densities or low curvatures. 
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Einstein's General Relativity (GR) describes the grav- 
ity field in a very successful way by the metric tensor of 
the space-time through the Einstein-Hilbert (EH) action. 
However, this action is non renormalizable and many au- 
thors have tried to solve different cosmological and as- 
trophysical puzzles by modifying it, i.e. with MGMs; for 
example, the introduction of Lagrange densities propor- 
tional to R n with n > 1 (Large MGMs or LMGMs) leads 
to Starobinsky inflation In the last years, terms pro- 
portional to R n with n < 1 (Small MGMs or SMGMs) 
have been taken into account to explain the present cos- 
mic acceleration without the typical cosmological con- 
stant (n — 0) or dark energy |2|. In early times, it is 
reasonable that terms with n > 1 can be significant but 
negligible today; whereas terms with n < 1 can be ne- 
glected in the early Universe but not in the present or 
late epochs, when the space-time curvature is very small. 

The introduction of new scalar curvature dependent 
terms in the action modifies Einstein's Equations (EEs) 
in the following way: 



[l + e'{R)]R llv --[R + e{R)]g llv 
+ 2 a p^V a ^ [e'(R)} = ^-, 

An) i 



(1) 



where l a piw = {g a pg^, v - g ai i9fiv) and e {n '(R) is the n 
derivative of e(R) with respect to the curvature. We are 
supposing a small modification of the EH action: 



SG =2 



d 4 ., 



-gM^ [R + e(R)} 



(2) 



with the reduced Planck mass Mpi = (SttGn)^ 1 ^ 2 
2.4 x 10 18 GeV (we are adopting the convention: Ti 
c = 1). We are focusing on SMGMs of the form: 



e(R) = A/i 2 



l R ri 



A = ±1. 



n < 1 (n ^ 0) . (3) 



However, other terms present a similar behavior (as a 
logarithmic one) and some results of this work are general 
for any e(R). To the above gravitational action, we have 
to add the standard matter one: 



Sm = 



M • 



(4) 



In the case that there is no matter, due to the symmetries 
of the vacuum, the solutions are maximally symmetric 
space-times with a constant curvature scalar i?o, which 
is solution of the following equation: 



R [1 - e' ] + 2e = , 



(5) 



where = e^(Ra). It is interesting to note that with- 
out a cosmological constant, Rq = 0, is not in general a 
solution for SMGMs. In fact, a large number of works 
have studied such modifications with fj, ~ Hn~ 10 -33 eV 
to explain the observed cosmic acceleration y, ll> IM III > 
whose origin could come from different compactifications 
of extra dimensions in M-theory • 

On the other hand, the third term in the left hand side 
of Q can produce important instabilities |8|, whose so- 
lution seems to demand the introduction of higher order 
terms (n > 1) in the scalar curvature 0]. However, in the 
present literature, it is possible to find three other prob- 
lems of SMGMs which are directly related to our present 
discussion [j| 0, . In this work we show that these 
three problems are connected and that their analyses are 
not enough to exclude SMGMs. On the contrary, we ar- 
gue that this type of models presents a well defined NL 
at intermediate scales. 

The NL in SMGMs was studied in 11] through weak 
field expansions around maximally symmetric vacuum 
solutions. It means that we can perform a series in the 
s(R) function supposing analyticity in the background 
solution R : e(R) = e + s' Q (6R) + e'^{5R) 2 /2 + ... with 
SR = R— Rq. For instance, taking trace in we obtain 
a covariant expansion, whose zero order is given by @ 
and the first order by: 



[1 + e' - Roe%] (SR) - 34'V Q V Q (8R) = - 



T 



(6) 



Here, T = g^T^ is the standard energy momentum 
tensor associated to (@J . On the other hand, if we perform 
a small perturbation of the background metric g^ v — 
9o ftu + 2$ (5° <5° , we can write the linearized fourth order 
equation for the metric perturbation $, and interpret 
this perturbation as the usual Newtonian potential. In 
|llj |. it was shown that we can recover an approximate 
Poisson equation and the NL if e' ' is small enough, i.e. at 



2 



distances r > (£q) 1/2 - But in the SMGMs under study: 

e"(R) = \n(n - 1){R/ fj, 2 ) 71 ' 1 /R. (7) 

Models with n < 2 (n ^ 0, 1) have typically large values 
of £q due to their inverse dependence on the Hubble scale 

On the other hand, in [10|, the gravitational force 
which appears in the paradigmatic model (n = — 1) 
due to a diffuse source in a locally de Sitter background 
was calculated. A linearly growing behavior was found, 
which is unacceptable because, for instance, it can in- 
crease the interaction between the Milky Way and An- 
dromeda by six orders of magnitude |l(Tj . To arrive at 
this conclusion, the authors also performed a similar ex- 
pansion and analyzed the first order tensorial equation 
from JU whose trace is given by ©. 

Finally, in Reference ||, it was shown that these 
SMGMs are equivalent to Scalar- Tensor Theories (STTs) 
excluded by Solar System (SS) experiments. The MGM 
action leads to fourth order equations, which can be stud- 
ied with the usual EEs if we add a new Scalar Degree of 
Freedom (SDF)-It is the so called Einstein frame fl2T|-. 
This new field has the following potential: 



M^[R(cp)e' v 



2[1 



(8) 



where R(f) is the solution of the equation e'(R) = 
exp{y/2/3ip/Mpi) - 1, and e£> = E^(R(<p)). Its mass is 
related to the second derivative of the potential: 



d 2 V v 1 R(<p)[3 - e' v ] + 4e ¥ 
dip 2 



3[l + e'J 2 



(9) 



In the vacuum solution, this field satisfies R(po) = Rq ~ 
fi 2 , which implies that the typical mass is of order m v ~ 
fx. However, if m v ~ Hq < 10~ 18 eV, the model can 
be excluded by SS tests, such as the deflection of light 
by the Sun, because ip mediates a new force with a long 
range dEa. 

In conclusion, several authors have detected important 
problems to the viability of SMGMs. Indeed, these three 
works are related because they have obtained inappro- 
priate behaviors of the gravitational theory taking into 
account its vacuum state. 

Different solutions have been proposed for some of 
these three problems For instance, Dick ^lj has pro- 
posed a fine tuning to save the gravitational potential. 
A model that verifies e'q = has the correct NL, but 
even more significantly, this tuning resolves the two other 
problems: the usual gravitational interaction between 
galaxies is recovered, at least at the linearized equations; 
and the SDF has a divergent mass ©■ Assuredly, the 
scalar field is not well defined on the vacuum state. This 
result can be understood because this field takes into ac- 
count the new degree of freedom of the gravity due to its 



fourth order equations, but if e'q = 0, the metric has as- 
sociated the standard second order equations in vacuum. 

A fine tuning usually means a physical misunderstand- 
ing, but there could be a fundamental reason for this can- 
cellation. For instance, Dick has proposed the following 
example \Ti\ : 



s(R) = -15(i 4 /R + 25n 6 /R 2 



(10) 



The vacuum solution reads Rq = 5/x 2 . These two terms 
have the same importance as the EH one in the vac- 
uum state, which means that the gravitational coupling 
is modified. The effective Mpi can be deduced from 
©: M 2 S — Mp,(l + e'q), and in this particular model: 
M| ff = 5M| 1 /6E3. 

On the other hand, in Reference [1J|, it has been 
shown that we can recover Newton's Gravity Law (NGL) 
through the Schwarzschild solution inside a de Sitter (or 
anti-de Sitter) space. In fact, the metric: 



ds 2 
A(r) 



-A(r)dt 2 
^ 2m 



AiryUr 2 
Rpr 2 
12 



l dn 2 



(11) 



is solution of the equations of motion QJ for a point- 
like source of mass m and with R given by Equation 
0. Therefore, we can argue that with a small enough 
background curvature any astrophysical test of gravity, 
which depends on the Schwarzschild solution and its NL, 
will be unaffected by the studied SMGMs. 

However, this solution is not completely satisfactory 
because it seems in contradiction with the previous re- 
sults. Furthermore, without a well defined NL, we can 
not identify the Mp\ in the action. In the same Ref- 
erence ^ nas been commented that in the ^ — > 
limit, we should recover the standard NL. Assuredly, if 
we remove the modification in the action, i.e. s(R) — > 0, 
we recover the EH one. However, if we follow the argu- 
ments presented in we arrive at the conclusion that 
the effective Mpi is given by 5Mpy/6 in the model 11U|) 
even in the fi — * limit. We can find similar surprises in 
the discussed results of 0,0]. Without the tuning, the 
three problems present more inappropriate behaviors in 
the EH limit. 

In our opinion, the question is that SMGMs modify 
the gravity at low energies. Therefore, it is consistent 
that the gravitational interaction was very different at 
low curvatures or low densities. This is the expected 
behavior, and the three discussed results agree with it 
showing that the gravitational interaction is more differ- 
ent for smaller values of \i. The interesting question is 
in what kind of physical environments we can use their 
results to contrast with experiments or observations. 

As it has been pointed out in 15], the expansions in 
the scalar curvature perturbations around the de Sitter 
vacuum solutions for the analyzed SMGMs JSJ present 
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the following form: 



e(R) = e 



(12) 



The experiments which provide us with the most precise 
values of NGL and the Mpi are realized on the Earth 
with SR/R ~ 0(pAir/#o M pi) ~ 1027 (taking into 
account the air density pAh- ~ 10 15 eV 4 ), which means 
that we can not identify the Mpi with the formula given 
in Reference |llj . The expansion can not be truncated 
inside the SS either; for example, we can roughly esti- 
mate the minimal SS curvature by the typical local Dark 
Matter (DM) value p^ 1 - 1CT 1 GeV/cm 3 - 1(T 6 eV 4 
\rj\. which implies 5R/R ~ 10 6 . 

Indeed, we can use this number as a typical value inside 
a galaxy. In such a case, we can not estimate the force 
between galaxies with the result reported in ^(J- I n fact, 
as their authors recognized, the linearized approximation 
around the vacuum solution inside the source (for exam- 
ple, the Milky Way) is not legitimate. Furthermore, the 
situation is analogous inside the body which can feel this 
possible force (such as Andromeda). Therefore, the cor- 
rect metric solution of these non linear equations, which 
drives its dynamic, can differ very much from the result 
found in 



> cry: 

m 



Finally, in the equivalent STT, we can not use the vac- 
uum values for the scalar field to study its effects in the SS 
(lif. Without performing any calculations, we can guess 
that SMGMs could present a very different gravitational 
interaction close to the vacuum. The Reference ja| has 
shown it in a very elegant way, but the experiments or 
observations which are used to constrain STTs are not re- 
alized on vacuum. The SS tests restrict severely STTs if 
the force mediated by the new scalar field has associated 
a long range. However, as in the chameleon case |18U19]. 
the characteristics of this field depend significantly on 
the local curvature and matter content, to which it is 
strongly coupled. The situation is again tricky, and to 
illustrate it, we can estimate the range of this new force 
inside the Sun supposing a Yukawa-type potential expo- 
nentially suppressed by m v . If R 3> /i 2 , we find by using 
Equations © and JJJ): 



m 



d 2 V v 
dip 2 



XR 



3< 



3n(n — 1) \fi 



(13) 



10 12 eV 



For instance, for the case n — —1: m v 
(SR/R ~ 10 30 with p Sun ~ 10 18 eV 4 ), we observe that 
the range is of order of m^ 1 ~ 10~ 19 m and conclude that 
the scalar field can not produce a force with observable 
effects out of the Sun (with a radius of r$ ~ 10 9 m). The 
situation is similar inside the Earth and all the typical 



SS gravitational sources 1 . 

The above estimations show that it is difficult to expect 
that the truncated expansions around the vacuum state 
were able to describe any physical observation. Practi- 
cally, this type of expansions can only be linearized for 
cosmological studies and only in very recent (or future) 
times when SR/Rq ~ Qm/^a ~ 0.3 |2(| (or smaller), but 
this is one of the most interesting features of SMGMs. 
They can modify the usual gravity at low curvatures 
and explain the present cosmology without a dark en- 
ergy component. 

If these models make any sense, they have to recover 
NGL at intermediate energies. We can not perform an 
expansion around the vacuum state, but we can perform 
an iterative method to resolve the Equation 0J. The 
zero order should be given by EEs (e(R) — ► 0): 



-Re fj,u — -^R^g'E^v 



T 



(14) 



Therefore, we find the usual relations: -Re = — T/Mpj 
and i?E M „ = (T^ - T g Etlv /2)/M^. To obtain the first 
correction to EEs, we can interpret the new terms like a 



new source: e 
we can write: 



(n) 



An)(_ T / M 2j 



Up to this first order, 



Ri 



(15) 



with 



T\ uu — 



M 2 l£E 



Te' 



pi 



5E fiv Z-eXixv 



(16) 



Here, the metric is given by EEs (|14fl and2Ea/3/ii/ = 
(5E Q /33E/^ - gE af j.gE/3u)- This procedure is inappro- 
priate for conformal matter (T=0) inside the studied 
SMGMs ©, but it is very helpful for non relativistic 
matter without pressure, T^„ = p<5 9 <5 9 : 



Mp> 2 



l + n 



(17) 



Clearly, we do not find any modification to EEs for 
e(R) = 0. Indeed, for a small enough value of /i, the 
correction is always negligible (n < 1). It is interesting 
to estimate when the maximum value of the new energy 
contribution (||Xi M „|| = MaxlTi^^l) is much smaller than 
p, which implies that NGL is recovered. For instance, if 



1 It is interesting to note that we find a negative square-mass term 
for the model with A = n = — 1 Q reproducing the instabilities 
studied in Q from an independent approach. 
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we can neglect the temporal and spatial variations of the 
source, we find the following relation: 

II^VIK IIT^II =► P ^fi 2 M^ : (n<l). (18) 

We observe that the same estimations, which were used 
to show that the vacuum expansions were incorrectly 
truncated, allow us to conclude that we have a well de- 
fined NL inside the SS for the proposed models with 
/i ~ H (n < 1). In fact, the two approximations are 
opposite in certain sense. This NL works for large densi- 
ties (or curvatures) in relation to the typical cosmological 
constant scale: M^Hg ~ 1CT 12 eV 4 (H§ ~ lCT 66 eV 2 ), 
whereas the linearized vacuum expansions work for small 
densities (or curvatures) with respect to the same scale. 

It is fair to say that we are not able to assure NGL for 
the interesting physical environment studied in . The 
Condition ljTS)l can not be satisfied in intergalactic spaces 
for fi ~ Hq. Indeed, we think that this is another very 
interesting property of SMGMs. They modify NGL at 
low densities, which is the expected behavior for n < 1. 
For example, if we neglect the mentioned temporal and 
spatial variations of the source, we obtain the following 
correction for the model with n = — 1 Q: 

a 4 M 4 

Tl ^ = \!±J^%%. (19) 

From this point of view, we have a new source of gravi- 
tation, which is certainly a DM source since we can only 
detect an anomalous metric behavior. It may be interest- 
ing to try to explain galactic dynamic anomalies or even 
rotation curves without DM halos, which requires the in- 
troduction of new unobserved particles in the standard 
framework [21). Similar approaches to general MGMs 
can be found in previous works with promising results 
0. However, if ||Ti M „|| > HX^H, higher order correc- 
tions have to be taken into account and a numerical cal- 
culation seems necessary. In this sense, we propose an it- 
erative method by repeating the described procedure up 
to an eventual convergence. Alternatively, if we know the 
metric, we can infer the necessary source. If this source 
does not agree with the observed one, we can deduce a 
possible function e(R) that can explain the difference. 
For simplicity, in a first step, we can study the functions 
given by Equation J3J) parameterized by fx, n and A. 

Before concluding this article, we would like to com- 
ment briefly about SMGMs inside the Palatini formalism, 
in which the metric and the connection are taken as inde- 
pendent variables We think that our discussion can 
also clarify several aspects concerning this approach. For 
instance, it is possible to find in the literature different 
works with different results about the NL in this formu- 
lation |l5ll22j|. In our opinion, it is fundamental to work 
at intermediate energies, as it has been performed very 
recently in (ljj . and not in the traditional low density 
approximation. 



We have studied the viability of SMGMs, which could 
be interesting to explain not only the present acceleration 
of the Universe but also its DM content. The idea is that 
they deviate from GR at low energies or low curvatures. 
In contrast, we recover the NL at intermediate scales. 
More analyses to clarify the possible existence of other 
kinds of problems have to be performed before postulat- 
ing them as a serious alternative. Work is in progress in 
this direction. 
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